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ABSTRACT 



A one-dimensional mixed layer model (Garwood, 1976, 1977) is used to 
simulate the mixed layer depth and temperature observed during POLEX, a 
component of the North Pacific Experiment. POLEX occurred during Janu- 
ary-February 1974 with mid-ocean observations (35'’N, 155‘’W) of tempera- 
ture and salinity made from R/P FLIP, while FLIP was under free-dri ft 
conditions. The results of this simulation show that apparent horizon- 
tal advection, due to the drift of R/P FLIP, was important in the heat 
and salinity budgets, but that vertical mixing and mixed layer depth 
changes were controlled primarily by the one-dimensional response of the 
turbulent kinetic energy budget to the local atmospheric forcing. Oc- 
casionally, surface salinity flux, due to large precipitation rates, can 
significantly alter (on the order of tens of meters) the thickness of 
the ocean surface turbulent boundary layer, with a demonstrated decrease 
in thickness during a single period of strong precipitation. 
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Buoyancy 

2 -3 

Turbulent buoyancy flux (cm sec ) 

Degrees Centigrade 

Non-dimensional drag coefficient 

Heat capacity of sea water (cal gm"^C~^) 

-2 -1 

Rate of evaporation (gm cm sec ) 

Vapor pressure of air 

Saturated vapor pressure of marine air 

Coriolis parameter 

Acceleration due to gravity 

Turbulent boundary layer depth (modeled) 

Mixed layer depth (modeled) 

Observed mixed layer depth 
Meters 

-2 -1 

Rate of precipitation (gm cm sec ) 
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Temperature in the air 

Sea-surface temperature 

Turbulent temperature flux (mwatts cm ) 

Time (hours) 

Frictional velocity at ocean surface 
Mean wind speed (m sec**^ ) 

Vertical turbulent velocity flux 
Velocity vector (cm sec~^) 

Velocity vector due to ocean currents 
Velocity vector due to Ekman flow 
Velocity of current due to quasi-geostropic flow 
Velocity due to the drift of R/P FLIP 
Vertical component of turbulent velocity 
Depth (meters) 

Depth increment 

Expansion coefficient for heat 
Expansion coefficient for salt 
Density of sea water (gm cm ) 

.O 

Density of air (gm cm ) 

Surface stress (dynes cm ) 

Temperature 

Represent the total instantaneous value 
Denotes the vertical average 
Increment of change in a variable 
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I. INTRODUCTION 



A. PURPOSE OF STUDY 

The objective of this research was to investigate the conclusions 
which were reached by Simpson and Paulson [1979] from ocean data ob- 
tained as part of the Pole Experiment (POLEX), a component of the North 
Pacific Experiment (NORPAX). POLEX was part of the first process ori- 
ented NORPAX experiment. The experiment was named POLE to indicate that 
the horizontal extent of the sampling was small compared to the largest 
scales investigated in NORPAX. Simpson and Paulson concluded: "One- 

dimensional mixed-layer deepening models failed to predict the mixed- 
layer depth and temperatures observed during POLE. Horizontal advec- 
tion, as evidenced from the salinity maximum frequently occurring at the 
bottom of the mixed layer and other near-surface changes in salinity and 
temperature not associated with local surface forcing, are responsible 
for the failure".^ 

This study will investigate the one-dimensionality of the POLE 
temperature and salinity data and determine if the mixed layer depth 
changes in response to the local atmospheric conditions (one-dimensional) 
or if these changes are caused by two and/or three dimensional advective 
processes. To accomplish this task, a non-stationary , one-dimensional 
bulk model of the mixed layer originally proposed by Garwood [1976, 1977] 

^ Simpson, J. J. and Paulson, C.A., 1979: Observation of Upper Ocean 

Temperature and Salinity Structure During POLE Experiment, J. Phys. 
Oceanoqr. , 9, page 869. 
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is expanded by including a salinity budget. A comparison of the re- 
sults of a strictly one-dimensional model of the upper ocean with obser- 
vations can be used to make an assessment of the one-dimensionality of 
the different variables which are involved in the upper ocean boundary 
layer response to the local atmospheric forcing. 

One hypothesis to be examined is that horizontal advection may be 
important for the heat and salt budgets, and at the same time vertical 
mixing and mixed layer depth may be controlled primarily by one- 
dimensional or local forcing alone. One objective of this study will 
be to determine how responsive are the variables of heat and salt in 
the upper ocean boundary layer to one-dimensional effects. A second 
objective of this study is to use the data obtained as part of POLEX to 
see if the one-dimensional mixed layer model proposed by Garwood will 
predict the observed mixed layer depth, even if advection dominates the 
local changes in heat and salt content in the upper ocean. 

The data used in this report contain observations made from R/P FLIP 
while under free-drift conditions. As a result of the free-drift condi- 
tion of FLIP during the data collection period, advection as seen by R/P 
FLIP will be a combination of the actual advection in an Eulerian coor- 
dinate system for the area and advection due to the drift of FLIP. 
Henceforth, the term advection will be used to mean a combination of the 
actual advection term and apparent horizontal advection due to the drift 
of FLIP. 

The data were collected during the period 28 January 74 through 
14 February 74. During that time, FLIP occupied a station approxi- 
mately 800 miles north of the Hawaiian Island Chain under free drift 
conditions. The position of FLIP ranged from 35“’39'N, 155‘’05'W to 
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34‘’36'N, 155*’25'W. The experimental site is hydrodynamically complex, 

as shown in Figure 1. The Subtropical Front is known to meander between 

32*’ and 35*’N (Roden, 1974; Barnett, 1976) . The region of the trade winds 

northeast of Hawaii has air-sea fluxes of latent heat in excess of 35 
-2 -1 

cal cm day (Wyrtki, 1965). The Subtropical water mass formed in this 
region contrasts markedly with the less saline Eastern North Pacific 
Central Water characteristical ly encountered north of SS^N. The Horse 
Latitudes are located only 3* of latitude to the south of the observa- 
tion area and the North Pacific Current is expected to affect the gen- 
eral hydrography of the region. 

B. MILITARY IMPORTANCE 

Mixed layer depth is an important factor when dealing with sound 
propagation in the ocean. It is used to determine such sonar detection 
variables as convergence zone distances, source level depth, and cutoff 
frequency for surface duct propagation. A model which could predict the 
change in mixed layer depth in response to the local atmospheric forcing 
could be incorporated in the existing acoustic models, for example the 
Fast Asymptotic Coherent Transmission (FACT) and the Naval Interim Sur- 
face Ship Model (NISSIM), presently being used at Fleet Numerical Ocean- 
ography Center (FNOC). A one dimensional mixed layer model would help 
the naval tactician to directly correlate the variability of the ocean 
with hisweapon system and determine the most efficient deployment of 
the naval Anti-Submarine Warfare (ASW) assets. 

C. LITERATURE REVIEW 

The vertical fluxes of heat, salt, and momentum across the sea-air 
interface are the sources of almost all oceanic motions. In the fully 
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turbulent oceanic mixed layer bounded by the sea-air interface above and 
by the dynamically stable water mass below, the vertical fluxes are 
large. Below the mixed layer, they are usually negligibly small so that 
one can decouple the mixed layer from the underlying stable, quiescent 
water mass. This homogeneity is the root of the term "slab", which is 
often used to describe the layer as depicted in Figures 2 and 3. 

When dealing with sound propagation in the ocean, the term mixed 
layer depth is sometimes taken to be that area with a uniform tempera- 
ture structure, but this is not always the case. In Figure 4 the 
bottom of the mixed layer depth (Region I) appears to be at about 80 
meters, and below there is a strong temperature gradient. However, ex- 
amination of the profile shows that the mixed layer depth is about 30 
meters. This shallow halocline was probably caused by an earlier rain- 
fall which decreased the near-surface salinity. The assumption of ver- 
tical homogeneity in a bulk model, because of the problem associated 
with solving for the interior fluxes of buoyancy and momentum, reduces 
to the need to know only the surface and entrainment fluxes. However, 
only small vertical gradients in these mean variables may be associated 
with large turbulent fluxes. Therefore, the slab assumption should not 
be as automatically applied to the turbulent kinetic energy budget. 

The immediate local reaction of this mixed layer to those fluxes re- 
sults in a homogeneous water column, i.e., vertical uniformity in the 
mean velocity and density. This adjustment to the density and velocity 
structure of the surface layers of the ocean to variable fluxes has been 
the subject of a large number of studies since Ekman's [1905] treatise 
where he originated the concept of a depth of frictional resistance for 
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the upper section of a wind stressed ocean. He further suggested that 
the wind-driven current will have a similar depth, the Ekman depth. 

Much of the one-dimensional theory for the ocean surface turbulent 
boundary layer or mixed layer seems to be dependent upon the validity 
of two basic hypotheses. The first of these is that vertical mixing 
within the turbulent boundary layer and entrainment mixing at its base 
occur in response to local atmospheric forcing, i.e., the surface_wind 
stress and net heat radiated or fluxed across the sea-air interface. The 
second hypothesis is that the mechanical energy budget is the key to the 
understanding and prediction of the mixed layer dynamics. 

Most observations show that local atmospheric forcing such as an 
increase in wind speed increases the mixed layer depth and low wind speed 
and high solar radiation decreases mixed layer depth. This is less true 
with increased depth below the mixed layer. Observations by Pollard and 
Tarbel [l975] showed that horizontal coherence between the fluctuating 
wind-driven currents near the surface was large and consistent with the 
local amplitude and phase differences of the wind history. With in- 
creased depth, the horizontal coherence decreased monotonical ly as did 
the vertical coherence of the deeper current measurements with those 
near the surface. Kroll [l975] pointed out that the time scale of these 
strongly inertial currents is on the order of a couple of days, depend- 
ing upon the stability of the underlying water column as well as the 
difference between the driving frequency and the local inertial fre- 
quency. The turbulent motions within the mixed layer are expected to be 
even more responsive to the local atmospheric conditions because of the 
small dissipation 'time scale for the turbulent kinetic energy--minutes 
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to hours, which depends upon the layer depth and turbulence intensity. 
Both the mean and turbulent kinetic energies in the mixed layer are ex- 
pected to reflect the local boundary conditions for wind stress and sur- 
face buoyancy flux. Measurements of turbulent velocity fluctuations, 
which were made by Powell et al [l975] and Jones and Kenney [1977], show 
that the root-mean-square turbulent velocity does scale with the water 
surface friction velocity, u*. Camp and Elsberry [l978] and Elsberry 
and Camp [l978] demonstrated the role of strong synoptic scale wind 
events for deepening. Garwood and Hal pern [l976] and Elsberry and Raney 
[1978] indicated the importance of diurnal heating, especially at times 
between storms. Hindcasting simulations by Thompson [1976], Mel lor and 
Durbin [l975] and Camp [l976] support the importance of surface boundary 
fluxes. 

Kraus and Turner |l967j improved the one-dimensional model originally 
developed by Kraus and Rooth [l96l]. This was accomplished by consider- 
ing the turbulent kinetic energy budget utilizing the heat equation and 
a mechanical energy equation. Because the frictional generation of heat 
is negligible, the vertically integrated heat equation provides a rela- 
tionship for the conservation of potential energy. They parameterized 
the mechanical production rate in terms of the friction velocity, but 
neglected the viscous rate of dissipation and the effects of entrainment 
shear production. 

The application of mechanical energy principles to the problem of 
mixed layer dynamics has caused some disagreement. An important contro- 
versy in mixed layer modeling has been centered around whether mean kin- 
etic energy or turbulent kinetic energy leads to the downward entrainment 
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buoyancy flux at the base of the turbulent boundary layer. The mean 
kinetic energy advocates, including Thompson jj97^ , Pollard 097?3 and 
Pri ce [1977J , argue that this mean flow energy is directly converted to 
potential energy by a mean flow instability that causes the bulk 
Richardson number to be less than some critical constant of the order 
one. Turbulent kinetic energy advocates do not disagree with the pos- 
sibility of this mechanism, but they argue it is not the dominant mech- 
anism because turbulence-generated instabilities will usually suffi- 
ciently erode the interface to preclude the occurrence of the mean flow 
instability. De Szoeke and Rhines jj975j and Yun 09783 demonstrated 
that wind-shear production of turbulent energy is ultimately the most 
important source of energy for deepening the mixed layer, if the surface 
buoyancy flux is neglected. Even when the mean flow energy is important, 
they find that the rate of conversion of this energy to potential energy 
via entrainment shear production is controlled by the constraints of the 
turbulent kinetic energy budget. The Kraus-Turner turbulent energy model 
has not proven to be entirely satisfactory, and a number of changes have 
been suggested. Stevenson 09793 examined a number of alternative para- 
meterizations for dissipation. Zeman and Tennekes 097^, Elsberry et al 
jl97^ , Alexander and Kim 09763 depth-dependent 

dissipation terms that helped to reduce unrealistic deepening* 

Garwood jj976, 1970 developed an ocean mixed-layer model using the 
Navier-Stokes equation of motion with the geostrophic component eliminated, 
the continuity equation in incompressible water, the heat equation from 
the first law of thermodynamics, the conservation of salt equation, and 
a linearized equation of state. The entrainment hypothesis depends upon 
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the relative distribution of turbulent energy between horizontal and 
vertical components and is offered as a plausible mechanism for govern- 
ing both entrainment and layer retreat. 

Garwood suggested that planetary rotation influences the dissipation 
of turbulence for deeper mixed layers and enables a cyclical steady state 
to exist on an annual basis. Furthermore, the rate of entrainment for the 
stable regime can not be a simple linear extrapolation of the unstable 
situations. Unlike the atomospheric boundary case, most of the solar 
radiation does not penetrate the layer. Therefore, downward turbulent 
heat flux in the oceanic boundary layer is as important as the upward 
flux during the course of both diurnal and annual cycles. The non- 
linearity of the interface entrainment tendency parameter, which is 
greatest for stable surface boundary conditions, results in a modulation 
of the long-term trend of the mixed-layer depth by the diurnal component 
of surface heat flux. In this model, buoyant production is somewhat 
more efficient than shear production as a source of energy for vertical 
mixing because of its unique effect on the vertical component of the 
turbulent velocity. The buoyancy equation is generated from the heat 
and salt equations together with an equation of state. 



In equations 1 and 2 0 is temperature, S salinity and p density 
while a and B are the expansion coefficients for heat and salt, respec- 
tively, and g is gravity. The tilde represents the total instantaneous 




( 1 ) 



and the definition for buoyancy. 




( 2 ) 



16 



value and the subscript zero denotes a representative but arbitrary con- 



stant value. The generalization of using b rather than 0 will cast the 
model equations in a form equally applicable to those situations where 
evaporation and precipitation contribute significantly to the surface 
buoyancy flux and the structure of the evolving pycnocline. The buoyancy 
equation also has a more obvious and direct role in the mechanical energy 
budget as shown in Figure 5. 

Miller 0976^] included salinity in the mixed-layer model of Kraus and 
Turner [l967] . He concluded that salinity is important in determining 
the density structure when applying a mixed-layer model to the entire 
ocean, but the relative importance of salinity effect on the short term 
evolution of the density profile may not be significant. Paulus 097^ 
included salinity in the mixed layer model of Camp [l97^. He concluded 
that the effects of salinity structure were noticeable only in deepening 
regimes. 
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II. DATA ANALYSIS 



A. METHODS 

The upper ocean temperature and salinity data and the meterological 
data during the NORPAX POLE experiment were obtained from Simpson and 
Paulson [ 1977 ] . The parameters required for surface boundary condition 
computations, which are used in Garwood's model, include wind speed, wind 
direction, cloud cover, sea surface temperature, air temperature (dry 
bulb), dew point and rain code. These parameters were recorded at ir- 
regular time intervals. A linear interpolation technique was used to 
obtain values for these input parameters at regular time intervals. The 
resulting interpolated data were compared with the original data. Fig- 
ures 6 through 12 show that no noticeable error was introduced into the 
parameters by the use of a linear interpolation. 

Temperature and salinity profiles as in Figures 3 and 4 were digi- 
tized by using a Houston Instrument Hi Pad Digitizer and a Tektronix 
4052 computer. The digitized profiles of temperature and salinity were 
linearly interpolated into hourly intervals. This provided for both 
initial conditions in the one-dimensional model and for fields to be 
used to compute observed mixed layer depth, h^(t), and to compare with 
subsequent model computations of T(z,t) and S(z,t). 

B. SURFACE FLUXES 

The mixed layer model uses the total surface heat flux (Q-j-(O)), in- 
cident solar radiation, wind speed at a height d, rate of evaporation 
(E), and rate of precipitation (P) to calculate surface fluxes of 
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buoyancy (heat and salt) and momentum. The friction velocity, in air, 



is calculated using the following formulas: 






u* = (T 5 /P 3 ) 



1/2 



(3) 

(4) 



where u^ is the mean wind speet (m/sec), 

a 



-3, 



c^ is the non-dimensional drag coefficient (1.4x10 ), 

-3 3 

p, is the density of air (1.25x10 gm/cm ), and 

a 

2 

T^ -is the surface stress (dynes/cm ). 

The turbulent fluxes of latent heat (Q^) and sensible heat (Q^^) are 
estimated using the following bulk aerodynamic formulas: 

(5) 






Qv, = - TJu 



a' a 



( 6 ) 



The net back radiation (Qj^) is estimated from the following empirical 
formula reported by Husby and Seckel [ 1975 ]: 

Q. = 1.14 X 10"^(273.16 + T (.39 - .OSeI^^) (1. - .60^) (7) 

D S u 

where is the saturated vapor pressure of the marine air directly in 
contact with the sea surface (.98 corrects for salt defects), 

T. is the air temperature (degrees Centigrade), 

a 

E, is the vapor pressure of air at approximately 10 m based on dew 

a 

point temperature, 

T^ is the sea-surface temperature (degrees Centigrade), and 
C is the fractional cloud cover. 

The upward heat flux is then: 



Q = Q + Qu + Qk 

^a ^e ^h ^b 



(8) 
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( 9 ) 



The solar insolation, Q^, was estimated by: 

Qs = (1. - ao:^') (1. - .66 C^)Qq 

The constants a and b are adapted from Tabata j_196^ and the cubic 
cloud cover correction from Laevast . The coefficient C is the 

fractional cloud cover and is the mid-day elevation angle of the sun. 

The clear sky radiation, Qq, is given by 

Qq=Aq+Ai cos(J)+Bi si n(|>+A 2 COs 2 (j>+B 2 S i n2(j) (10) 

from the formula developed by Seckel and Beauday [19733 The coeffi- 
cients (Aq, etc.) were calculated by harmonic representation of the values 
presented in the Smithsonian Meteorological Tables (List, 1958) with 
(j)=(2'n-/365) (t-21 ) where t is the Julian day of the year. Turbulent 
temperature flux, 

01 ) 

} pCp 

turbulent velocity flux, 

^(0) = 'ii (12) 

turbulent salinity flux, 

SV(0) = {P-E)S,„, (13) 

and turbulent buoyancy flux, 

b^(0) = g [cT^(0) - es^(o)] (I'D 

where a is the coefficient of thermal expansion, 

3 is the coefficient of contraction of sea water due to salinity 
variations, and 

is the surface salinity value. 
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Figures 13, 14, 15 and 16 are graphs which show the hourly values of 



wind stress, x, total heat flux at the surface, Qj, salinity flux, S'w', 
and buoyancy flux, b'w' , respectively, as calculated by the one-dimensional 
model. These figures can be compared with Figures 17 and 18 which show the 
interpolated hourly values of wind stress, x, total heat flux at the sur- 
face, Qj, buoyancy flux, b'w' , and various other components which were re- 
ported by Simpson and Paulson ^197^. The values of these fluxes cal- 
culated here are nearly identical to the values reported by Simpson and 
Paulson. In Simpson and Paulson's graph on buoyancy flux (Figure 17), a 
positive buoyancy flux results when > 0. In the one-dimensional 

model (Figure 16) a positive buoyancy flux results when < 0 and E > P. 

1 . Surface Wind Stress 

The wind stress calculated using the bulk approximation with a 

-3 -2 

drag coefficient of 1.4x10 has a mean value of about 0.66 dynes cm 

for the experimental period. Figure 13 shows the output for the surface 

wind stress as calculated by the one-dimensional model. Analysis of 

Figure 13 shows that for periods up to as long as 4 days the wind stress 

_2 

is considerably less than 1 dyne cm , and it is evident that only two 
periods of relatively high winds were encountered during the experiment. 

One period of high wind near the middle of the experiment (time = 200-260 

-2 

hours) with a maximum wind stress of 2.5 dynes cm , and the other at the 

_2 

end of the experiment with a maximum wind stress of 3.3 dynes cm . 

2. Surface Heat Flux 

Comparison of the total heat flux with net all -wave flux (Figures 
17 and 19) shows that the net all-wave flux is generally the dominant 
term in the surface heat balance. However, a two-day period of enhanced 
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I 

I 




latent heat flux, due largely to a significant increase in wind speed 
midway through the experiment (time = 250 hours), is evident when com- 
paring Figures 13 and 14. 

Simpson's [ 1977 ] calculation of the total heat budget for the 

_2 

experiment suggested that the ocean gained 201 cal cm for the time per- 
iod from 3 Feb 74 - 12 Feb 74. Simpson's heat budgets calculations indi- 
cated that ocean gained heat for 6 days, was in near thermal equilibrium 
with the atmosphere for 2 days and lost heat to the atmosphere for 2-day 
period. These results, coupled with the low wind stress, imply that the 
upper dynamics might have been dominated by net surface heating for a 
significant part of the experiment. 

3. Salinity Flux 

Figure 15 shows the salinity term as computed by the one- 
dimensional model (equation 13) without including the precipitation term 
in the computation. The precipitation term was omitted in the calculations 
because no reasonable correlation could be made between a rain code given 
in the experiment data report (Simpson and Paulson 1977) and an amount of 
precipitation. The results of model runs which included precipitation 
(estimated by the rain code) were compared to model runs which neglected 
precipitation altogether. These model comparison runs showed no signifi- 
cant difference except during one period when precipitation was much 
larger than evaporation. Simpson and Paulson jj979^ observed that over- 
all evaporation and precipitation were in near equilibrium, with precip- 

_2 

itation exceeding evaporation by 23 mg cm . 

4. Buoyancy Flux 

Stratification in the ocean exerts a strong dynamical influence 
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on mixing since a stably stratified fluid requires work against buoyancy 
forces if mixing is to occur. For a stably stratified fluid, the density 
increases with depth and the buoyancy flux is directed downward, while 
the converse is true for an unstable, stratified fluid. In the one- 
dimensional model, surface buoyancy flux, b'w‘ , is determined by the 
air-sea transfers (equation 14). Figure 16 shows the model computed 
buoyancy flux with < 0 and E > P resulting in a positive buoyancy 
flux. The net heat flux. Figure 14, dominates the surface values of 
buoyancy flux during the experiment. During daytime periods, negative 
buoyancy flux results from the domination of the solar heating at the 
surface, while during the night periods a positive buoyancy flux is 
produced by the combined effect of net long-wave radiation and the up- 
ward turbulent fluxes of heat and moisture. 



C. ADVECTION 



The determination of the importance of advection to the heat and 
salt budget is necessary before an assessment can be made on how one- 
dimensional is the response of the variables of heat and salt are in 
the upper boundary layer. The velocity vector V of the advection term 
(V'(VT,VS)), is separated into two components in order to determine 
which component dominates the advection term. 

1 . Heat and Salinity Budgets 

Neglecting internal sources and sinks, the equation for conser- 
vation of heat or salt is: 

0 



a 



3(T,S)+ rv(T,S) + 3(T,S)'w' + 3(T,S)'u' + 3(T,S)'v' 
3t 3z 3x 3y 

(a) (b) (c) (d) (e) 



dz=0 (15) 
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where T is heat, S is salinity, t is time, d is depth and V is the velo- 
city vector in a coordinate system moving with the drifting platform. 

The meanings of these terms are: 

(a) Heat or salinity change with time, 

(b) Heat or salinity change due to advection in the moving coordinate 
system, assuming that when averaged over a day such features such 
as internal waves would not contribute to advection, 

(c) Vertical heat or saltflux divergence, and 

(d) and (e) Horizontal heat or salinity flux divergence. 

Terms (d) and (e) are assumed to be small because they can not be 
evaluated. Therefore, equation (15) can be rewritten in terms of the 
variable temperature: 

dz = - (T^o)^ ^ ^(d) 

(f) (g) 

Term (f) is equal to the surface heat flux Qj/pCp, where Qj is the net 

-2 -1 

surface heat (Qj“Qs“*^b“*^e“*^h^*^c^ ’ density of sea 

-3 -1 -1 

water, 1 gm cm , and Cp is the specific heat of sea water, 1 cal gm C 

Term (g) is assumed to be negligible if d is well below the mixed layer, 

i.e., d > The vertical integrated mean advection of temperature is 

therefore, 

■/ 

-d 
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The same procedure used to derive (17) may be followed to find the 



average advective change in salinity. The surface salinity flux is 
given in equation (13). Thus the vertically integrated mean advection 
of salinity is: 



A 



VS)dz = 



-d 



r 0 - 

r 




r 0 n 
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Sdz 
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Sdz 


J 

*--d 
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/d . 



At 






S(E-P)dt 



(18) 



The magnitude of actual advection (in an Eulerian coordinate system) 
can be checked by estimating V and VT or VS. The term V(VT, VS), is the 
advection quantity averaged over a time interval. At. The velocity vector 
V can be separated into two components, 

V = V^H-Vp (19) 

where is the velocity vector due to the ocean currents and is the 
velocity vector due to the drift of R/P FLIP. The velocity vector due to 
the ocean currents is itself composed of two components: 

V, = V ^ (20) 

c g e 

where Vg is the velocity of the current due to the quasi -geostrophic flow 
and is assumed to be parallel to the surface temperature contour. There- 
fore Vg'VT is assumed to be negligible. The order of magnitude of V^, 
the velocity current due to Ekman flow, is a function of the wind stress 
and can be estimated by, 

IVgl = IT^I = 1.57 cm sec‘^ = 0.0565 km hour'^ (21) 

"W 

where is the average wind stress, 0.66 dyn cm"^, p is the density of 
sea water, h is average mixed layer depth, 5000 cm and f is the Coriolis 
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parameter, 2cosiri(j), where co is the angular velocity of the earth and (}> 
is the latitude. 

The velocity calculated due to Ekman flow is insignificant when com- 
pared with the drift rate velocity calculated for R/P FLIP (Vp=0.3 km hour"^) 
from 28 Jan 1974 to 14 Feb 1974. The calculated Ekman velocity is not suf- 
ficient enough to explain the drift of POLE. Hence, R/P FLIP was by no 
means flowing with the same mixed layer water mass throughout the experi- 
ment. Therefore, the velocity vector V will depend on the drift rate of 
R/P FLIP, V = Vp. What Simpson and Paulson |j979^ call advection was not 
actual advection for the area but apparent advection due to the drift of 
R/P FLIP. 

2. Calculation of the Rates of Advection 

The distance R/P FLIP traveled from 0700 hours 28 Jan 1974 to 
0700 hours 14 Feb 1974 (408 hours), was calculated to be 121 km, using 
a rhumbline distance calculation and a position of 155°05'W to 

34‘’36N, 155°25’W. Vertical profiles of temperature and salinity were 
taken from R/P FLIP from 0000 hours 30 Jan 1974 through 0700 hours Feb 
1974 (total of 367 hours) which covered a distance of approximately 
109 km. Figures 19 and 20 show the observed changes in heat and salt 
fluxes. The amounts of advection in a coordinate system moving hori- 
zontally with R/P FLIP were then estimated by subtracting the accumu- 
lated surface heat and salt fluxes from the changes in heat and salt 
content, respectively. Figures 21 and 22 show the results of this cal- 
culation. Least square fits were performed on the results, giving a 

-2 -1 

slope of 2.0 cal cm hour for heat content gain and a slope of 
0.45 *^/00 cm hour"^ for salinity gain. These values can be used to 
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compute the average horizontal temperature and salinity gradients by 
multiplying by the total hours of observation (367 hours) and dividing 
by the average mixed layer depth (approximately 40m) and the distance 
FLIP traveled during the observational period (109 km), 

|V^T1 = 0.0017 C km‘^ (22) 

|V^S| = 0.0004 °/00 km'^ (23) 

where |V^T| is the horizontal temperature gradient and |V^S| is the hori- 
zontal salinity gradient. 

For comparison purposes, gradients were obtained from a NORPAC 
Atlas jj96^ by dividing the atlas temperature change (AT = 0.2 “’C) and 
the atlas salinity change (AS = 0.05 ^/OO) between FLIP starting and 
final positions (35“’39'N, 155°05'W to 34'’36'M, 155‘’25'W) by the total dis- 
tance (121 km) which FLIP traveled, 

|V^T|^ = 0.0017°C km‘^ 

IVp^Tl^ = 0.0004 °/00 km'^ 

where and |V^S|^ are the horizontal temperature and salinity 

gradients that were calculated from the NORPAC Atlas. Although these 
values agree with the earlier calculation, equations (22) and (23), they 
were only obtained to show agreement in order of magnitude with recorded 
atlas values. These results further illustrate the fact that what 
Simpson and Paulson [j97^ call advection was probably the result of the 
drift of R/P FLIP. 

Figures 19 and 21 show that the temperature field, which is a 
function of both time and depth, has large fluctuations which can only 
be explained by advection, or the drift of FLIP which is a non one- 
dimensional effect. These large fluctuations dominate the heat budget 
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on the time scales shorter than 100 hours, but over a longer time scale, 
surface flux is equally important (one-dimensional effects). Figure 23 
is a plot of observed sea surface temperature and a plot of the model - 
computed mixed layer temperature, with the model -computed temperature 
slightly offset on the temperature scale for comparison. Comparing the 
two plots, a one-dimensional component can be detected in the observed 
sea surface temperature time series but the stronger component due to 
the drift of FLIP is evident. Therefore, this apparent advection (as 
seen by FLIP) of temperature is important on all time scales, but domi- 
nates the short time scale. Figures 20 and 22 show that advection of 
salinity dominates the salinity budget over both the short and long time 
scales up to at least 400 hours, which was the period of observations. 

D. MIXED LAYER DEPTH 

1 . Observed Mixed Layer Depth 

The observed mixed layer depth, shown in Figure 17 as used by 
Simpson and Paulson jjOToJ, is the shallowest depth at which the den- 
sity is not more than 0.02 units greater than the density at a depth 
of 5 m. However, salinity and temperature changes can occur above 
z = 5m. In such cases, the above definition will result in a misleading 
computed time series for mixed layer depth. A more practical definition 
for observed mixed layer depth, which would be more consistent with the 
precision of the available POLEX observational data, would be the shal- 
lowest depth at which the observed value of is 0.02 units greater 
than the observed surface value. For comparison with the Garwood model, 
the observed mixed layer depth (hQ) will be defined according to this 
latter definition. 
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Contours of temperature for the observation period are shown in 
Figures 24 and 25. Figure 24, which is based on a contour interval of 
0.5 C, shows a mixed layer depth as varying from near the surface to ap- 
proximately 90 meters. Figure 25, which has a contour interval of 0.05 C 

2 

after having been smoothed twice, shows a mixed layer from near the sur- 
face to approximately 80 meters. This illustrates the importance of 
resolution when determining mixed layer depth from temperature contours. 

A higher resolution frequently results in a shallower mixed layer depth. 
Therefore, a temperature contour interval of 0.5 C is not sufficient 
resolution to accurately determine mixed layer depth from the contour 
plot. Although density is a function of both temperature and salinity, 
during this time series the salinity influence was not significant, so 
that the density was mainly dependent upon temperature. A mixed layer 
depth definition could then be defined as the shallowest depth at which 
the temperature is not more than 0.1 C less than the temperature at the 
surface which would compare well with the similar definition based on the 
observed density structure, as is evident by comparing Figure 25 with the 
observed mixed layer depth in Figure 26. 

2. Comparison of the Model -Computed Mixed Layer Depth with Observed 
Mixed Layer Depth 

Figure 26 is a time series of the depth of the model computed 
turbulent boundary layer, -h(t), and of the observed mixed layer depth, 
-hQ(t). The value of h should be less than or equal to the mixed layer 
depth as determined from observed 3p/9z if the model simulation of the 

2 

Smoothing was accomplished by giving the center data point a weight 
of 2 and all ,of the surrounding data points a weight of 1. 
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vertical density structure changes are accurate. Figure 27 is a plot of 
the same functions, h(t) and hQ(t) between time = 200 hours and time = 400 
hours. Notice that the model results compare favorably with the observa- 
tions except for a period of about one day starting at hour 285. Figure 
28 is the model computed mixed layer depth, h-j(t). The value of -h-j(t) 
is taken to be the depth at which the model is 0.02a^ units greater than 
the value of <a^> where < > denotes the vertical average for the turbu- 
lent boundary layer. This definition is equivalent to that for hQ(t) if 
the turbulent boundary layer is homogeneous. Figure 29 shows both h.j(t) 
and the observed mixed layer depth, h^Ct) for the same time period as 
Figure 28. Notice that h-|_< h for all t. Although hQ(t) compares well 
with h(t), there are times during which hQ(t) is more comparable to h.|(t). 
For example, for the period between t = 285 hours and t = 310 hours, 
hg(t) is most similar to h-|(t). This happens because, although the tur- 
bulent boundary layer may be quite shallow, the associated transient 
thermocline is not sufficiently strong to influence the values for h-j(t) 

and hg(t). 

Only during one time period, at about time = 100, is the observed 
mixed layer depth found to be significantly deeper than either the tur- 
bulent boundary layer, h(t), or the mixed layer depth, h-j(t), computed 
by the model. This difference between the predicted and observed mixed 
layer depth could be the result of R/P FLIP drifting through a warm core 
eddy or filament of water mass having more tropical characteristics. 
Figures 21 and 22, which are the graphs of the advective components of 
the heat and salt budgets, show a large increase in the magnitudes of 
both the advection of heat and salt at about time = 100 hours. This 
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large increase in advection of both heat and salt is consistent with the 
hypothesis that R/P FLIP drifted across either an eddy or a filament of 
water from the south of the Subtropical front. The Subtropical water 
mass which lies to the south is warmer and more saline than the Eastern 
North Pacific Central water. 

3. Salinity Effects on the Mixed Layer Depth 

An additional model simulation was conducted in order to deter- 
mine what effect precipitation has on the mixed layer depth. A steady 

-5 -2 -1 

rate of precipitation of 5 x 10 gm cm sec was prescribed over a 12- 
hour period starting at t = 200 hours. At all other times precipitation 
was assumed to be negligible. This event corresponds in magnitude, at 
least, with the rate of precipitation reported by Simpson [1977] for 
5 February 1974. Figure 30 shows a rapid decrease in the surface salin- 
ity with the injection of the precipitation. It drops from about 34.11 *^/00 
to 33.99 *^/00. This drop in surface salinity causes a rapid rise in the 
surface salinity flux as shown in Figure 31. Figure 32 shows the model 
prediction of the mixed layer depth with and without the precipitation 
event. The associated increase in downward buoyancy flux, b'w'^Qj due 
to the increase in salinity flux (Figure 31) caused the model mixed 
layer depth to deepen much less rapidly than for the case without pre- 
cipitation. The effect would have been even more apparent if the pre- 
cipitation period had not coincided with strong winds. 

Analysis of Figure 29 shows that the turbulent boundary layer 
depth predicted by the model was approximately 20 meters deeper than the 
observed mixed layer depth at t = 210 hours. This time period coincides 
with that reported by Simpson and Paulson [l979]| as the period of most 
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intense precipitation. These results suggest that, if the amount of 
precipitation were accurately known and included in the model boundary 
conditions, there would be an even better correlation between the model 
mixed layer depth and the observed mixed layer depth. 
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III. SUMMARY 



This analysis of 15 days of salinity, temperature and density pro- 
files from R/P FLIP as part of the NORPAX POLE experiment and the ac- 
companying model study suggest that, non one-dimensional effects were 
very important for both the heat and salt budgets. Vertical mixing and 
mixed layer depth are controlled primarily by the one-dimensional response 
of the turbulent kinetic energy budget to the local atmospheric forcing. 

It is apparent that the non one-dimensional effects may be largely ex- 
plained by the drift of R/P FLIP. Results from POLE site showed that this 
apparent horizontal temperature advection (as observed by the moving plat- 
form) was important on all time scales and dominated the small time scale. 
Nevertheless, the one-dimensional model computation of mixed layer depth 
compared favorably with the observed mixed layer depth. It is interest- 
ing that even though this apparent advection usually dominated the salt 
budget over both the long and short time scales, the one-dimensional mixed 
layer model response to the surface salinity flux can be significant. 
During the single period of strong precipitation, the mixed layer depth 
was shallower (on the order of tens of meters) than it would have been 
without precipitation. This is consistent with the results of Miller 
[l976) , showing that salinity flux can at least occasionally play an 
important part in mixed layer dynamics in the subtropics. Further 
studies would have to be conducted with observational data that included 
the time, duration, and amount of precipitation before a firm conclusion 
can be made on the longer-term importance of salinity effects on the 
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mixed layer depth. Qualitatively, model performance was improved when 
the influence of precipitation was taken into account. 

Simpson and Paulson [j979l concluded that the one-dimensional models 
failed to predict the mixed layer depth and temperature observed during 
POLE. Why did the models of Pollard et al j2973]| and Niiler ^197^ , which 
were used by Simpson and Paulson, fail to predict the observed mixed layer 
depth? One part of the explanation for the model's apparent failure could 
be the definition of mixed layer depth used by Simpson and Paulson. An 
accurate determination of the mixed layer depth from a vertical profile 
can only be made by considering the entire profile, an erroneous analysis 
of the mixed layer depth is probable whenever there is a transient 
pycnocline above z = 5 meters. Therefore, contrary to the results of 
Simpson and Paulson [1973 , it has been shown that the one-dimensional 
turbulence closure model used here does seem to explain most of the ob- 
served mixed layer depth changes and that the mixed layer depth is more 
influenced by the local surface fluxes than by horizontal advection. 
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Figure 1. The location of the R/P FLIP during the POLE 

experiment (35*N, 155‘’W) in relation to general 
oceanic features. (Simpson and Paulson 1977) 
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Figure 2. Idealized density and mean velocity profiles of the ocean 

mixed layer ( ). Mixed layer depth is (h), (6) is the 

thickness of the interface or entrainment zone. Region I 
is the fully turbulent mixed layer depth. Region II is the 
slightly stable, intermittently-turbulent entrainment zone. 
Region III is the stable underlying watermass having negli- 
gible vertical fluxes in comparison to those of Region I. 
(Garwood 1977) 
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Figure 3. Typical POLEX data temperature, salinity, and sigma-t pro- 
files taken from Simpson and Paulson [1979]. Added to the 
profile is Region I the ful ly-turbulent mixed layer of 
depth h. Region II the slightly stable, intermittently- 
turbulent entrainment zone of thickness 6. Region III is 
the stable underlying water mass having negligible fluxes 
in comparison to those of Region I. 



37 



SIGr:n-T 

I • 1 » \ i 1 1 ! * 1 

26.4 25.6 26.8 26.0 26.2 26.4 

SfiLIMITY in PPT 

I 1 1 1 i 1 1 1 1 

34.06 34.10 34.16 34.20 34.26 

TEMPERATURE IN DEO C 

j , 1 , J ; J , , , j 

10 11 12 13 14 16 




Figure 4. Typical POLEX data temperature, salinity, and sigma-t pro- 
files taken from Simpson and Paulson [1979]. Added to the 
profile is Region I the fully- turbulent mixed layer of 
depth h. Region II the slightly stable, intermittently- 
turbulent entrainment zone of thickness 6. Region III is 
the stable underlying water mass having negligible fluxes 
in comparison to those of Region I. 
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Figure 5. Mechanical energy budget for the ocean mixed layer. Asterisks 
indicate those processes that must be parameterized to close 
the system of equations. (Garwood 1977) 
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measured wind direction (solid line 
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Rain Code: 0-no rain, 1-fog, 2-mist, 3-light rain, 

4-rain. 
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Figure 13. Model -computed hourly values of the wind stress (t). 

The origin of the abscissa is 1000 GMT 28 January 1974, 
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Figure 14. Model computed hourly values of the total downward heat 
flux at the surface. The origin of the abcissa is 1000 
GMT 28 January 1974. 
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Figure 15. Model computed hourly values of salinity flux (S'w'), neglecting 
precipitation. 
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Figure 16. Model computed hourly values of buoyancy flux (b'w') 
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Figure 17. Interpolated hourly values of mass and buoyancy fluxes, 
surface salinity, sea surface temperature, surface 
sigma-t, and the various components of the heat flux 
involved in the determination of the buoyancy flux. 
(Simpson and Paulson 1977) 
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Figure 18. Interpolated hourly values of the various components 
of the surface momentum and heat balance are shown. 
(Simpson and Paulson 1977) 
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Figure 19. The observed changes in heat content in the upper 
40 meters (solid line). Model computed accumulated 
surface heat flux (dotted line). The origin of the 
abscissa is 0000 hour (local) 30 January 74. 
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Figure 20, The observed changes in salt content in upper 40 
meters (solid line). Model -computed accumulated 
surface salinity flux (T^), dashed line, which 
IS negligible in comparison with the observed 
changes in salt content. The origin of the 
abscissa is 0000 hour (local) 30 January 74. 
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Figure 21. Advection of heat (solid line) determined by 
subtracting the accumulated heat flux from the 
change in heat content (Figure 21), Least 
squares fit to data resulted in a slope 2.0 cal 
cm“^ hour" . 
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Figure 22. Advection of salt (solid line) determined by 
subtracting the accumulated surface salinity 
flux from the changes in salt content (Figure 
20). Least square fit to data resulted in a 
slope 0.45 0/00 cm hour"l. 
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Figure 23. Comparison of model computed mixed layer temperature (solid line) with 

observed for sea surface temperature (dotted line). Model computed mixed 
layer temperature has been offset slightly for the comparison. The origin 
of the abscissa is 1200 GMT 30 January 1974. 
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Figure 24. Isotherm distribution for the period 30 Jan 74 through 
14 Feb 74, The contour interval is O-S^C. 
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Figure 25. Isotherm distribution, smoothed twice, for period 30 Jan 74 
through 14 Feb 74. The contour interval is 0.05°C. 
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Figure 26. The dashed line with the symbol (0) is the observed mixed layer depth, 
hQ(t). The solid line is the model computed turbulent boundary layer 
depth, h(t). The origin of the abscissa is 1000 GMT 28 January 1974, 
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Figure 27. The dashed line with the symbol (0) is the observed mixed layer 
depth, hQ(t). The solid line is model computed turbulent 
boundary^layer depth, h(t). 
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Figure 28. Model computed mixed layer depth, h,(t). The origin of the abscissa 
is 1000 GMT 28 January 1974. 
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Figure 29. Model computed mixed layer depth (solid line), h-|(t). Observed ^ 
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Figure 31. Model computed hourly values of salinity flux (S'w‘ = (p-e)Sq) 
with hypothetical precipitation (5 x 10“^ gm cm"2 sec"l), 
starting at time = 220. 
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The dashed line shows the model computed mixed layer depth, h(t), 
without considering the effect of precipitation. The solid line 
shows a decrease in the model computed_^ixed la^er de|j»th as a result 

of including the precipitation (5 x 10 gm cm sec ). 
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